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We consider dimers on the star lattice (aka the 3-12, Fisher, expanded kagome or triangle- 
honeycomb lattice). We show that dimer coverings on this lattice have Z2 arrow and pseudo-spin 
representations analogous to those for the kagome lattice, and use these to construct an exactly 
solvable quantum dimer model (QDM) with a Rokhsar-Kivelson (RK) ground state. This QDM, 
first discussed by Moessner and Sondhi from a different point of view, is the star-lattice analogue of 
a kagome-lattice QDM analyzed by Misguich et al. We give a detailed analysis of various properties 
of the classical equal-weight dimer model on the star lattice, most of which are related to those of the 
RK state. Using both the arrow representation and the fermionic path integral formulation of the 
Pfaflian method, we calculate the number of dimer coverings, dimer occupation probabilities, and 
dimer, vison, and monomer correlation functions. We show that a dimer is uncorrelated with dimers 
further away than on neighboring triangles, that the vison-vison correlation function vanishes, and 
that the monomer- monomer correlation function equals 1/4 regardless of the monomer positions 
(with one exception). A key result in the fermionic approach is the vanishing of the two-point 
Green function beyond a short distance. These properties are very similar to those of dimers on 
the kagome lattice. We also discuss some generalizations involving arrow representations for dimer 
coverings on "general Fisher lattices" and their "reduced" lattices, with the kagome, squagome, and 
triangular kagome lattice being examples of the latter. 
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I. INTRODUCTION 

Classical dimer models (CDMs) have long been of in- 
terest in statistical and theoretical physics j^'^'^'^ An im- 
portant breakthrough was ushered in with the solution 
of the dimer problem on the square lattice in terms of 
Pfaffians.^'^'* In a subsequent seminal paper, ^ Kasteleyn 
generalized this solution method to the dimer problem on 
any planar lattice, and also showed explicitly that CDMs 
can exhibit singularities, corresponding to phase transi- 
tions, as a function of the dimer weights. A few years 
later Fisher showed that the Ising model on an arbitrary 
planar lattice can be mapped to a CDM on a different pla- 
nar lattice,!*^ thus making it possible to solve the former 
class of models with the Pfafhan method developed for 
the latter. More recently, Moessner and Sondhiii have 
demonstrated that the nature of the phases and phase 
transitions in the two models connected by this map- 
ping can be very different. Specifically, they revisited 
Fisher's mapping between the ferromagnetic Ising model 
on the square lattice and the CDM on the star lattice^ 
(shown in Fig. [T]) and found that the symmetry-breaking 
transition in the former model maps to a confinement- 
deconfinement transition for test monomers in the latter. 

Over the last two decades, quantum dimer models 
(QDMs) have come to play an important role in the 
search for phases with exotic types of orders and exci- 
tations in quantum antiferromagnets (for recent reviews, 
see Refs. [islandfl^. Most intriguingly, there is the possi- 
bility of spin liquid phases^^ with deconfined, fractional- 
ized (spin-1/2 spinon) excitations. These spin liquids do 
not break spin rotation or spatial symmetries but are in- 
stead characterized by unconventional quantum ordersr^ 



(e.g. topological order) not describable by local order 
parameters. The canonical example of a spin liquid is a 
short-range resonating-valence-bond (RVB) statei^ given 
by a macroscopic superposition of products of spin sin- 
glets formed by pairs of nearest-neighbor spins. QDMs 
were introduced by Rokhsar and Kivelson^- to be an ap- 
proximate, effective description of the low-energy physics 
when it is dominated by such short-range singlet fluctu- 
ations. 

The set of dimer coverings in a CDM is taken to form 
an orthogonal basis for the Hilbert space of QDMs on the 
same lattice. A QDM Hamiltonian can contain both po- 
tential and kinetic energy terms, which are respectively 
diagonal and off-diagonal in this dimer basis. If these 
terms have amplitudes v and t respectively, the QDM 
ground state will be a function of v/t. This parame- 
ter space generally includes a so-called Rokhsar-Kivelson 
(RK) point, at which the ground state is an equal- 
amplitude superposition of dimer coverings (in each topo- 
logical sector) As a consequence, there is a close rela- 
tionship between QDMs at the RK point and CDMs: 
any ground state correlation function at the RK point 
that only involves operators which are diagonal in the 
dimer basis (e.g., the dimer-dimer or vison-vison corre- 
lation function) is equal to the corresponding correlation 
function for the CDM with equal dimer weights on the 
same lattice. 

The search for dimer liquids in QDMs was crowned 
with success when Moessner and Sondhii^ showed that 
the ground state of the triangular-lattice QDM is a dimer 
liquid (i.e. has exponentially decaying dimer-dimer cor- 
relations) over a finite range of v/t. In another remark- 
able paper, Misguich et al}'^ introduced and analyzed a 
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QDM on the kagome lattice that differed in several re- 
spects from the QDMs studied until then. This kagome 
QDM does not have potential energy terms, and its ki- 
netic energy involves flippable loops^- of various lengths 
(in contrast, the prototypical QDM studied in Refs. 17 
and [isl only includes the shortest flippable loops). Due 
to special properties of the kagome lattice, which allow 
for Zi arrow and pseudo-spin representations of its dimer 
coverings Misguich el al)^ were able to show that 

this QDM can be written in an extremely simple form in 
terms of the pseudo-spin variables. As a consequence, 
all its eigenstates and eigenvalues are exactly known. In 
particular, its ground state in each topological sector is 
the RK (equal-amplitude) state. On the kagome lattice 
this state is an extremely disordered dimer liquid: the 
dimer-dimer correlations vanish for dimers further apart 
than on neighboring triangles. This has been shown from 
the arrow representation's, as well as from Pfaffian cal- 
culations of the dimer-dimer correlations of the CDM 
on the kagome lattice i ^^i^'* Another interesting result for 
the CDM on the kagome lattice is related to monomer 
deconfinement: the monomer-monomer correlation func- 
tion equals 1/4 regardless of the distance between the two 
monomers. Again, this has been shown both from the ar- 
row representation's, and from Pfaffian calculations.— 




FIG. 1: The star lattice. The lattice sites are the vertices 
in this figure. The lines connecting the sites will be called 
bonds. There are two types of bonds, which will be referred 
to as triangle bonds and linking bonds. A triangle bond is 
part of a triangle. A linking bond connects two neighboring 
triangles. 



In this paper we consider dimers on the star latticsi^ 
(Fig. [T|), a non-bipartite lattice that is one of the 11 
Archimedean tilings--^ Because of its frustration^^ and 
small coordination number (z = 3) this is an interest- 
ing lattice to consider in the search for unconventional 
phases in spin and dimer models. Various studies of this 
kind have recently appeared in the literature. We have 
already mentioned Ref. [Ill which we will come back to 
shortly. From an analysis of exact diagonalization spec- 
tra Richter et al^ concluded that the kagome and star 
lattices are the only Archimedean tilings for which the 
ground state of the nearest-neighbor spin- 1/2 Heisen- 



berg model does not have magnetic long-range order. For 
the star lattice the ground state is instead a unique va- 
lence bond crystal (VBC) with singlets across the bonds 
that link different triangles ;^^'^^ If the exchange across 
these linking bonds is made sufficiently small compared 
to that of the triangle bonds, exact diagonalization re- 
sults suggest that the ground state is a different, three- 
fold degenerate VBC.^° We also note that an experimen- 
tal realization of a star antiferromagnet, [Fe3(/i3-0)(/x- 
OAc)6-(H20)3][Fe3(Ai3-0)(Ai-OAc)7.5]2 -THaO, has been 
reported recently."^^ Furthermore, Yao and Kivelson^ 
have shown that the ground state of the Kitaev model 
on the star lattice is a chiral spin liquid. 

We show that dimers on the star lattice have proper- 
ties that are very similar to those for the kagome lattice 
summarized above. We first show that dimer coverings 
on the star lattice have arrow and pseudo-spin represen- 
tations analogous to those for the kagome lattice. We 
use these to construct a star-lattice QDM that has simi- 
lar properties as the kagome QDM discussed by Misguich 
et al^ In particular, its ground state is the RK state. 
This star-lattice QDM was first discussed in Ref. :l2 as 
a special case of a QDM introduced there via a modified 
Fisher construction mapping it to a three-dimensional 
Ising model in the time-continuum limit. We give a de- 
tailed analysis of the star-lattice CDM for equal-weight 
dimers (the case related to the RK ground state as men- 
tioned above), using both the Pfaffian method and the 
arrow representation to calculate the number of dimer 
coverings (on a finite lattice without boundaries), dimer 
occupation probabilities, and dimer, vison, and monomer 
correlation functions. 

For our Pfaffian calculations we use Samuel's formu- 
lation based on a path integral description of free Ma- 
jorana fcrmions.'^'^ We find that the (two-point) Green 
function of the fermions vanishes beyond a short distance. 
This is a key result that causes a dimer to be uncorre- 
lated with all dimers further away than on a neighbor- 
ing triangle. It furthermore makes the calculations of 
the vison- vison and monomer- monomer correlation func- 
tions tractable; although in the fermionic formulation 
these functions contain a "string" going from one vi- 
son/monomer to the other, which makes the number of 
pairings from Wick's theorem grow exponentially with 
the vison/monomer separation, the number of nonzero 
pairings is strongly reduced by the above property of the 
Green function. We also use the arrow representation 
to calculate the various quantities considered. These ar- 
row derivations are conceptually simpler and more intu- 
itive than the Pfaffian ones, in part because they lend 
themselves to a pictorial representation. Finally, we con- 
sider some generalizations to general Fisher lattices (any 
lattice hosting the dimers in Fisher's mappingiS), and 
discuss "reduced" lattices of these (examples of which 
include the kagome, squagome, and triangular kagome 
lattice). 

The paper is organized as follows. The arrow and 
pseudo-spin representation are introduced in Sees. [Til and 
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mil In Sec. mil we also construct the exactly solvable 
QDM. The number of dimer coverings is considered in 
Sec. HVl The Green function is calculated in Sec. |Vl 
In Sees. I VII and IVIII we discuss dimer occupation prob- 
abilities and dimer correlations. The vison-vison and 
monomer-monomer correlation functions are calculated 
in Sees. IVIIII and IIXI In Sec. |X] we summarize simi- 
larities and connections between dimers on the star and 
kagome lattice. In Sec. |Xl]we discuss properties of gen- 
eral Fisher lattices and their reduced lattices (with some 
technical details relegated to Appendix [B]) . Some con- 
cluding remarks are given in Sec. IXIII In Appendix |^ 
we present an alternative pseudo-spin representation of 
dimer coverings on the star lattice. 



II. ARROW REPRESENTATION OF DIMER 
COVERINGS 

In this section we show that there is a one-to-one cor- 
respondence between dimer coverings and arrow configu- 
rations on the star lattice. The arrows are Ising degrees 
of freedom located on the lattice sites. First, however, 
we establish some basic terminology that will be used 
throughout the paper. 

The star lattice consists of triangles and dodecagons 
as shown in Fig. [TJ The lattice sites sit at the vertices 
of the polygons. The thin lines connecting different sites 
will be referred to as bonds. While all sites in the star 
lattice have the same local environment (by definition of 
an Archimedean tiling), there are two types of incquiv- 
alent bonds, which will be referred to as triangle bonds 
and linking bonds. A triangle bond is part of a triangle 
(as well as a dodecagon) while a linking bond is part of 
two dodecagons and is a link between two triangles. 

A bond can be occupied by a dimer, which will be 
drawn as a thicker line on the bond. The dimer touches 
the two sites that are connected by the bond. We will 
consider close-packed, hard-core dimer configurations. 
Close-packed means that no site is left untouched by 
a dimer. Hard-core means that no site is touched by 
more than one dimer. The fact that each site is touched 
by exactly one dimer will be referred to as the closed- 
packed, hard-core constraint. In the following we will 
refer to such dimer configurations simply as dimer cov- 
erings. (In Sec. IIXI we will relax the close-packed part 
of the constraint to consider dimer-monomer configura- 
tions in which two sites in the lattice are not touched by 
a dimer.) An example of a dimer covering on the star 
lattice is shown in Fig. [2lItop). 

Elser and Zeng showed that dimer coverings on the 
kagome lattice allow for an arrow representation.^^ Mis- 
guich et al. made extensive use of this representation 
in their work on QDMs on the kagome lattice j^^i^^ and 
also pointed out that the arrow representation can be 
generalized further to all lattices made of corner-sharing 
triangles. We will now show that an arrow representation 
exists for dimer coverings on the star lattice too (arrow 



representations for some other lattices will be discussed 
in Sec. |Xl] and |XlTl) . 




FIG. 2; Top: Example of dimer covering on the star lattice. 
Dimers are shown as thick lines. Bottom: The same dimer 
covering in the arrow representation. 



In the arrow representation for the star lattice, an ar- 
row is located on each lattice site. It can point in one of 
two directions: either towards the center of the triangle 
to which the site belongs, or in the opposite direction, 
towards the center of the site's linking bond. In the first 
case the dimer that touches the site lies on one of the 
site's two triangle bonds. In the second case the dimer 
lies on the site's linking bond. If a triangle contains a 
dimer on bond ij the arrows on sites i and j both point 
into the triangle while the arrow on the triangle's third 
site points out of the triangle. In contrast, if a triangle 
has no dimer, all three arrows point out of the triangle. 
As these are the only possibilities, each triangle must 
satisfy the constraint that an even number (0 or 2) of 
its arrows must point into the triangle. Next consider a 
linking bond kl. If it has a dimer, the arrows on sites 
k and / both point "into" (i.e. towards the center of) 
the linking bond. Alternatively, the linking bond has no 
dimer, in which case the directions of both arrows are 
reversed to point "out of (i.e. away from the center of) 
the linking bond. Hence each linking bond must satisfy 
the constraint that an even number (0 or 2) of its ar- 
rows must point "into" the bond. The constraint on a 
triangle is exactly the same as for the lattices of corner- 
sharing triangles for which the arrow representation has 
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been used earlier. Furthermore, the new constraint in- 
troduced for Hnking bonds takes the same form as the 
triangle constraint: an even number of arrows must point 
inwards. An example of a dimer covering and its arrow 
representation is given in Fig. [2l 



III. PSEUDO-SPIN REPRESENTATION OF 
DIMER COVERINGS AND AN EXACTLY 
SOLVABLE QUANTUM DIMER MODEL 

Elser and Zeng introduced an Ising pseudo-spin rep- 
resentation for dimer coverings on the kagome latticci^ 
Misguich et alw^^ developed this representation further 
by unveiling additional important properties. In this sec- 
tion we show that an analogous pseudo-spin represen- 
tation exists for the star lattice. In Sees. IIII Al and 
IIIIBI we derive the properties of this representation. In 
Sec. IIII CI we use the pseudo-spin representation to de- 
fine a star lattice QDM whose eigenstates and energy 
spectrum are exactly known and which has a Rokhsar- 
Kivelson ground state. This model has previously been 
discussed by Moessner and Sondhi from a different point 
of viewiii 



A. The operators (t^[D) 

We consider a star lattice on a closed surface, i.e. a sys- 
tem without boundaries. The dimer coverings (DCs) can 
be grouped into topological sectors. To do this one looks 
at the transition graph of two DCs obtained by superpos- 
ing them on top of each other. These transition graphs 
may contain two different types of non-intersecting even- 
length closed loops: loops that enclose a nonzero area 
and loops that enclose zero area (loops of the latter type 
have length 2). Loops of the first type can be topologi- 
cally trivial or nontrivial (all loops of the second type are 
obviously topologically trivial). "^^ If the transition graph 
of two DCs has no topologically non-trivial loops, these 
DCs are in the same topological sector. 

Next consider a given topological sector for DCs on the 
star lattice. Let |co) be some arbitrary but specific DC in 
this topological sector; it will play the role of a reference 
DC. Introduce a spin-1/2 operator (T^(Z3) for each do- 
decagon D with eigenvalues c^iD) = ±1. To any DC |c) 
in the same topological sector as |co) we can associate a 
pair of pseudo-spin configurations (PSCs) {a^{D)} and 
{—a^{D)} related by a global spin reversal, as follows: 
Draw the transition graph (co|c). Consider the closed 
loops that enclose a nonzero area in this graph; these 
loops are taken to be domain walls separating dodecagons 
with opposite values of (T^.^° The need to identify the two 
PSCs related by a global spin reversal with the same DC 
is due to the fact that the domain walls only determine 
the value of the products a^{D)a^{D') for all dodecagons 
D, D' , which are invariant under a global spin reversal. 
The reference DC |co) is associated with the pair of PSCs 



having all up and all down. 

The mapping from DCs in a given topological sector 
to pairs of PSCs can also be formulated in terms of the 
arrow representation. Let D and D' be two neighboring 
dodecagons and consider the linking bond separating D 
and D'. If the two arrows associated with this linking 
bond point in opposite (the same) directions in |c) and 
I Co) , there will (will not) be a domain wall between D and 
D' in the transition graph (co|c), so that the product 
a''{D)a^{D') = -1 (1) in the pair of PSCs associated 
with |c). 

Two different DCs |c) and |c') in the same topological 
sector will map to different pairs of PSCs. To see this 
we note that if the DCs are different, there will be some 
linking bonds that are occupied by a dimer in |c) but 
unoccupied in |c') and vice versa. (In other words, it is 
impossible to have dimer coverings that differ only in the 
dimer occupation numbers on triangle bonds). Thus the 
arrows on these linking bonds point in opposite directions 
in |c) and |c'), and therefore the product a^{D)a^{D') for 
the dodecagons D and D' separated by each such linking 
bond will be different for the PSCs of |c) and |c'). 

Any pair of PSCs related by global spin reversal cor- 
responds to a unique DC |c) in the topological sector of 
I Co) (this DC will however depend on the choice of |co)). 
To see this we note that given these PSCs we can find 
the value of the product a^{D)a^{D') for each neighbor- 
ing pair of dodecagons D and D'. If the product is — 1 
(1) the directions of the two arrows on the linking bond 
separating D and D' are opposite (the same) in |c) and 
I Co). Since every arrow is associated with a unique link- 
ing bond in this way, all arrows are assigned a direction 
by this procedure, thus uniquely specifying |c). 

From this mapping between DCs in a given topologi- 
cal sector and pairs of PSCs {a^{D)} and {—(t^{D)}, it 
follows that the number of DCs in a topological sector 
of the star lattice is the same for each sector and given 
by 2^°~^ where No is the number of dodecagons. This 
result can be expressed in terms of the number of sites 
N and the genus g of the manifold by using Euler's for- 
mula V + F — E — X where V, F, and E are the number 
of vertices, faces, and edges in the lattice graph and the 
Euler characteristic x = 2 — 2g for a closed orientable 
surface of genus g. Here V = N, F — No + Nt where 
Nt = N/3 is the number of triangles, and E = Nt + Ni 
where Nt = 3Nt = N is the number of triangle bonds 
and Ni — N/2 is the number of linking bonds. This gives 
Nd — N/6 + 2 — 2g so the number of DCs in each topolog- 
ical sector is 2^/6+1-29. in Sec. HV^we will show that 
the total number of DCs on this lattice graph (i.e. includ- 
ing all topological sectors) is 2^^^+^. Thus the number 
of topological sectors is 2^/6+1/2^/6+1-29 ^ = 4^, 
which is the expected result for a non-bipartite lattice on 
a genus-g manifold. 
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B. The operators <t^(D) 

Using the arrow constraints, one finds that there are 
64 different configurations for the 18 arrows on the six 
triangles that surround a given dodecagon. Each of these 
arrow configurations corresponds, in the dimer picture, 
to a "flippable" dimer configuration along one of the 
32 even-length loops that enclose the dodecagon and an 
even number of its six surrounding triangles (for each 
such loop there are two flippable dimer configurations). 
There is 1 loop of length 12, 15 loops of length 14, 15 
loops of length 16, and 1 loop of length 18. The loops 
in these four classes surround the dodecagon and 0, 2, 4, 
and 6 of its neighboring triangles, respectively. 

We now define the operator (T^(Z?) as 

a-(i^)=^(|L„P))(L„P)|+h.c.). (1) 

a 

This definition is analogous to the one for the kagome 
latticCi ^^'^^ The sum runs over the 32 even-length loops 
associated with D. \La{D)) and \La{D)) are the two 
"flippable" dimer configurations along loop a around do- 
decagon D. It can be seen that, in the arrow represen- 
tation, (T^(Z?) flips the direction of the 12 arrows sitting 
on the dodecagon D. This operation respects the arrow 
constraints since it flips an even number (0 or 2) of the 
arrows involved in each arrow constraint in the system. 
From this result it immediately follows that ((T^(£)))^ = / 
(flipping the arrows twice is equivalent to doing nothing) 
and that [a^ (D), {D')] = (it doesn't matter in which 
order the arrows are flipped). 

Next let D' be any of the dodecagons neighboring 
D. Consider the effect of acting with a^{D) on a DC 
|c) in the same topological sector as |co). If the direc- 
tions of the two arrows on the linking bond separat- 
ing D and D' are opposite (the same) in |c) and |co), 
they will be the same (opposite) in (t^{D)\c) and |co). 
Thus for all D' that neighbor D, a^{D) anticommutes 
with a^{D)a^{D'). Furthermore, if Da and D}, are two 
neighboring dodecagons, both different from D, then 
<T^{Da)<7^{Db) clearly commutes with a^{D) since the 
latter doesn't affect the arrows on the linking bond sepa- 
rating Da and Dt,. These two results can be generalized^ 
to, respectively, D' not a neighbor of D and Da and Di, 
not neighbors of each other. From this we conclude that 
a^{D) anticommutes with (J^{D) and commutes with 
a^{D') for D' ^ D. This shows that (as implied by the 
notation) a^{D) is the spin-1/2 operator that flips the 
pseudo-spin <J^{D) at dodecagon D. Also note that 

\{a^{D)^I. (2) 

D 

This constraint involving all dodecagons follows from the 
fact that n_D ^'^{D) is the operator that effects a global 
spin reversal, under which a PSC {(T^(I?)} maps to its 
partner {— (t^(£))} that represents the same DC. Alterna- 
tively, this constraint can be understood from the arrow 



representation, since the Ihs of Eq. ([2|) flips all arrows 
twice and therefore is equivalent to doing nothing. 

Starting from the arrow representation for the refer- 
ence DC |co), and given one of the two PSCs for the DC 
|c) in the same sector, the arrow representation for |c) 
can be obtained either by applying to |co) at all do- 
decagons with o"^ = -|-1 in |c) or by applying o"^ to |co) 
at all dodecagons with — —1 in |c). Either way, the 
arrows that point in different directions in |co) and |c) 
will be flipped once, while the arrows that point in the 
same direction in |co) and |c) will be flipped zero or two 
times. 



C. An exactly solvable quantum dimer model 

Given the form of (t^(Z?), it is natural to consider a 
QDM defined by the Hamiltonian 

H = -TY,^^{D) (3) 

D 

with r a positive constant. This QDM, which has no 
potential energy terms, is the star-lattice analogue of the 
kagome lattice QDM introduced by Misguich et al}^ It 
has been discussed previously by Moessner and Sondhiii 
as a particular limit of a QDM they arrived at via a 
mapping (based on a modified Fisher construction) to 
a three-dimensional Ising model in the time-continuum 
limit. 

Let us discuss some basic properties of the QDM ([3]) 
from the point of view of the pseudo-spin representa- 
tion that we have introduced in this section. Clearly the 
eigenstates are of the form Y[d W^i^)) where \a^{D)) 
is an eigenstate of a^{D) with eigenvalue a^{D) — ±1. 
Note however that due to the constraint only eigen- 
states with an even number of dodecagons in the excited 
CT^ = — 1 state are allowed. These cr^ = — 1 excitations 
localized at the dodecagons are (point) visons.^^'^^ We 
will discuss them further in Sec. IVIIII In the ground 
state, each dodecagon is in the = +1 state which is 
just the sum of the cr^ — +1 and = —1 eigenstates. 
Expanding out the product over dodecagons, the ground 
state is seen to be the equal-amplitude superposition of 
all dimer coverings in a given topological sector, i.e. the 
RK state. Since the ground state energy is the same in 
each of the 4^ topological sectors the system has topolog- 
ical order li^i^ (Note that the ground state is degenerate 
even for a system with a finite number of sites.) 

The fact that the ground state is 43-fold degenerate 
can also be seen as follows: A ground state |5'o) of ([3]) 
satisfies o'^(Z?)|4'o) = |^o) for each dodecagon D. This 
can be regarded as a constraint on |5'o)- There is one 
such constraint for each dodecagon. However, due to ^ 
only N]j — 1 of these constraints are independent. Every 
independent dodecagon constraint reduces the number of 
possible states in the ground state manifold by a factor 
of two, since it specifies one of two possible eigenstates 
at that dodecagon. As in Sec. IIII Al we invoke the result 
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to be shown in Sec. IIV Al that the total number of dimer 
coverings is 2^/^+^ which is therefore also the dimension 
of the Hilbert space. Thus the dimension of the ground 
state manifold is 2^/^+-'^~(''^''~^) = 4^, where we used 
Nd = + 2 - 2g as shown in Sec. HITH 

As noted in the introduction, RK-state correlations in- 
volving operators that are diagonal in the dimer basis are 
given by the corresponding correlation functions in the 
CDM with equal dimer weights. This CDM will be ana- 
lyzed in the following sections. 



IV. THE NUMBER OF DIMER COVERINGS 

In this section we consider the number of dimer cover- 
ings -2 on a star lattice graph with a finite number of sites 
N . We will take the graph to be embedded on a closed 
surface so that the system has no boundaries. We first 
consider a surface of arbitrary genus g and give a deriva- 
tion of Z based on the arrow representation introduced 
in Sec. |TT1 Next we calculate Z for a torus (genus g = 1) 
using the Pfaffian method. Some alternative derivations 
of this result are discussed in Sec. IXII and Appendix 1X1 



A. Arrow derivation for a genus-g manifold 

Consider a star lattice graph embedded on a closed 
orientable surface of genus g. If this graph has N sites 
there are triangles and iV/2 linking bonds. As the 
two latter numbers must both be integers, A'^ is a mul- 
tiple of 6 (and thus also even). Since there is one arrow 
per site and each arrow is an Ising variable (i.e. it can 
point in two different directions) it follows that in the 
absence of any constraints there would be 2^ different 
ways to choose the directions of the arrows. Each of the 
triangles and linking bonds does however contribute an 
arrow constraint, giving a total of A^/S -I- A^/2 = 5A^/6 
constraints. As will be shown below, one of these con- 
straints can be deduced from the others, leaving 5A^/6— 1 
independent constraints, each of which can be used to 
eliminate exactly one arrow variable."*^ Thus the total 
number of dimer coverings on this star lattice graph is 

Z = 2^-(5^/^"i) = 2^/^+1 (4) 

We note that the result is independent of the genus g. 
The associated entropy S = log Z is thus 

S^^^N + \og2. (5) 
6 

As shown in Sec. IIII Al there are 4^ topological sectors, 
each of which have the same number of dimer coverings 
Zts = Z/A^. The corresponding entropy per sector is 
therefore 

5ts = ^^+(l-2g)log2. (6) 



From these expressions one sees that the entropy per 
site in the thermodynamic limit is given by (1/6) log 2, 
in agreement with previous calculations.— In addition 
to the leading 0{N) (i.e., extensive) term the entropy 
also has a sub-leading 0(1) term that carries informa- 
tion about the topology of the system. 

Let us show that only 5N/ 6 — 1 arrow constraints are 
independent as claimed above. Consider an arbitrary site 
i. Denote the triangle and the linking bond it belongs to 
by T and I respectively. Define the arrow variable a^^T = 
+1 (—1) if the arrow on site i points out of (into) T. 
Similarly, define the arrow variable Ui^i = +1 (—1) if the 
arrow points out of (into) I. Clearly these two variables 
are related by ai^T = —Uij. Furthermore define ax = 
riiGT "^^.T where the product runs over the three sites 
belonging to the triangle T, and a; = Yiiei '^i^ where the 
product runs over the two sites belonging to the linking 
bond I. We have 

(n (n = n = (-i)^ = i- (7) 

On the other hand, the 5A^/6 arrow constraints can be 
written 

aT^ai = l for all T, I. (8) 

From ([7]) we see that one of the constraints in ([8]) fol- 
lows from the others, leaving bN/6 — 1 independent con- 
straints. 



B. Pfaffian calculation for a torus 

The Pfaffian method can be used to calculate the num- 
ber of dimer coverings Z on any planar lattice graphj^ 
In this method the central object is an antisymmetric 
matrix A, sometimes called a Kasteleyn matrix, whose 
dimension equals the (even) number of sites N. The 
matrix element Aij is nonzero only if sites i and j are 
connected by a bond in the lattice graph, in which case 
the magnitude of Aij is set equal to l4^ The sign of Aij 
defines the direction of an arrow^ on the bond ij: If 
Aij > the arrow points from i to j. The signs of the 
matrix elements of A should be chosen to satisfy Kaste- 
leyn's clockwise-odd rule:^-''* The perimeter of each face 
of the lattice graph should contain an odd number of ar- 
rows pointing in the clockwise direction around the face. 
One then has Z = |Pf A|, where Pf A is the Pfaffian of 
A. Using (Pf Ay = det A it follows that Z = VdetA for 
a planar lattice graph. 

Fisher^° and Wu- have used the Pfaffian method to 
calculate the free energy per site/dimer for the dimer 
problem on an infinite star lattice. For our Pfaffian cal- 
culations we follow Fisher in viewing the star lattice as 
a square lattice of 6-site unit cells, as shown in Fig. [3l 
where we also use his labeling of sites and choice of bond 
arrow directions to satisfy Kasteleyn's sign rule. A site 
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X. X. i.. 

-ix-1,y;^^: x,y :^^;x+1,y;^^ 



FIG. 3: The star lattice drawn as a square lattice (left) of 
6-site unit cells (right). Also shown are the directions of the 
Kasteleyn arrows and the labeling of unit cells and sites within 
a unit cell. These conventions follow Ref. 



i is written i — r,a where a = 1, . . . , 6 is the site label 
within the unit cell, and r = (x, y) labels the unit cell, 
where x — l,...,Nx is the horizontal coordinate and 
y = 1, . . . , Ny the vertical one, with and Ny being 
the number of unit cells in the horizontal and vertical 
direction, respectively. The total number of sites in the 
lattice is then N = 6NxNy, which is always even so that 
the lattice can be completely covered by dimers for any 
choice of Nx and Ny . 

When the sites on opposite edges in the two directions 
are connected by additional bonds to form a torus, the 
resulting lattice graph is no longer planar. However, as 
shown by Kasteleyn;^ the dimer problem on a torus can 
still be solved using Pfaffians by a small modification of 
the procedure described so far. In this case Z can be 
expressed as a linear combination of the Pfaffians of four 
different Kasteleyn matrices A^'^'^y {v^, Vy — 0, 1), where 
the first (second) superscript on A^''^^ labels the type of 
boundary condition imposed on the matrix elements in 
the horizontal (vertical) direction. Compared to the sin- 
gle Kasteleyn matrix for the open-boundary case, these 
four matrices differ only in the presence of additional 
nonzero matrix elements corresponding to the additional 
bonds introduced to establish the torus geometry. Specif- 
ically, the new matrix elements of A'^^"" are given by 



with A: 



"ljV,,a,5;l,y,l 



(_1)^. 



for all y, 
for all X, 



(9) 
(10) 



.^^ j . From the arrows in Fig. [3] one 
sees that the value (1) for Vx, Vy implies periodic (anti- 
periodic) boundary conditions on the matrix elements. 
On the torus one then has^ (see also Ref. ,49.) 



(11) 



where r^^,^^ are signs whose values will be determined 
below. 

In order to calculate these determinants, as well as 
other quantities to be considered later, we will use 



Samuel's reformulation of the Pfaffian method in terms 
of a path integral description of noninteracting Majorana 
fermions living on the lattice sites.— A fermion on site i 
is associated with a Grassmann variable tpi . The action S 
of the fermions is defined in terms of a Kasteleyn matrix 
A of the dimer problem. 



(12) 



The partition function Z of the fermions then equals the 
Pfaffian of A, 



Jv^ exp(-S') = Pf A. 



(13) 



Here / 'Dip = J dipi ■ ■ ■ dipN- In the following the matrix 
A is taken to be one of the four Kasteleyn matrices A''^"" 
for the dimer problem on the torus. To avoid cluttering 
the notation we will however not write the superscripts 
explicitly. We first introduce new Grassmann variables 
Tpk^a related to the Tpr,a by a Fourier transformation. 



y k 



Ak-r 



(14) 



The form of the wavevectors k depends on the bound- 
ary conditions Vx and Vy, i.e. kx — 27r(n^ -I- Vx/2)/Nx 
and similarly for ky. The integers Ux and Uy are cho- 
sen so that the sum over k runs over the first Bril- 
louin zone. The partition function can now be written 
Z = jVij) exp(-S') with 



(15) 



k a/3 



where A{k) is the 6x6 matrix 



/ 





1 


1 












-1 





1 














-1 


-1 





1 


















-1 





1 


1 














-1 





1 




V 










-1 


-1 





J 



(16) 



From it then follows that dot A = Z^ is given by 



det A = JJdeti(fc). 



(17) 



Because det A{k) = 4 for all k, det A = 4^-^" = 4^/*^ in- 
dependently of the boundary conditions Vx and Vy . Thus 
(fTTjl becomes Z — (r/2)-2^/^ where r = roo + ?'oi + ?'io + 
rii. We expect r to contribute to an A'^-independent term 
in the entropy (cf. Eqs. ©-([H])). Thus r can be deter- 
mined by manually counting the number of dimer cover- 
ings for a torus with only one unit cell"'^ (i.e. N = 6). 
We find Z = 4 for this case, which implies r — A, i.e. all 
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terms in (jlip come with a positive sign. Thus the num- 
ber of dimer coverings on a star lattice on a torus with 
iV sites is given by^ 

Z = 2^/6+1, (18) 

in agreement with the resuh ([4]) derived from the arrow 
representation in Sec. IIV Al 

V. THE GREEN FUNCTION 

Within the framework of the fermionic formulation of 
the Pfaffian method;^ quantities of interest for the dimer 
problem can be expressed in terms of fermionic correla- 
tion functions.-^iSiiSaiSO The average number of dimers on 
a bond is given directly as the two-point Green function 
of the Majorana fermions. Furthermore, dimer, vison, 
and monomer correlation functions can be expressed in 
terms of multi-point Green functions, which in turn, by 
Wick's theorem, are given as sums of products of two- 
point Green functions. In this section we calculate the 
two-point Green function of the star-lattice CDM for the 
case of equal dimer weights. 

The expectation value of an arbitrary operator O in 
the fermionic system defined by and is (O) = 



P = 



Q = 



R = 



Since the dependence on r, r' in Eq. (pO| is only 
through r — r' = R, we define Gr,a-r',a' = G{R)aa' — 
G{R^,Ry)aa'- We see from Ecis. ' ([2T])-(I111) that each 
term in the matrix elements of A~^{k) contains at most 
one power of e^"^^^ and at most one power of e^*'^" . From 
Eq. (HOI) it then follows that 

G{R^,Ry)a,a' ^0 if > 2 or if > 2, (25) 

i.e. the Green function vanishes identically beyond a very 
short distance. The same property has recently also been 
found for the Green function on the kagome latticei^^i^i 
For each R, G{R) = G{Rx,Ry) defines a 6 x 6 ma- 
trix. Of these, due to Eq. only G(0,0), G(1,0), 
G(-1,0), G(0,1), G(0,-1), G(-l,l), G(l,-1), G(l,l) 
and G(— 1, —1) can be nonzero. However, not all of these 



Z ^ J Vij] O exp(— 5). The two-point Green function 
(henceforth just called the Green function) is then 

G,j = (V-.Vj) - (19) 

Thus, viewed as a matrix, the Green function is the 
inverse of the Kasteleyn matrix: G = A^^ . Using 

Eq. (|14p . one can write Gr,a-r',a' = (V'r.aV'r'.Q') = 

[N.Nyy^j:,^ ,^, -^'^iik^Jk',..'). FromEq. m 

one finds {i)k,a4'k' ,a') = Sk'~kA~^{k)aa', where A^'^{k) 
is the inverse matrix of A{k). Thus 

G..„;w.a' ^ r ^ r ^e^'^<^-^')A-\kU, (20) 

where we have also taken the thermodynamic limit Nx, 
Ny —^oo. In this limit the Green function becomes in- 
dependent of the boundary conditions i^x , i^y because the 
wavevectors are no longer discrete. The inverse matrix 
of A{k) takes the form 

where P, Q, and R are 3x3 matrices given by 



(22) 



(23) 



(24) 

I 

matrices are independent because G = A^^ is antisym- 
metric and thus G{—Rx,—Ry) = —G{Rx,Ry)^- 

The results for the Green function obtained in this sec- 
tion will be used in the fermionic calculations of the quan- 
tities considered in Sees. IVIIIIXl 



VI. DIMER OCCUPATION PROBABILITIES 

Let i and j be two sites connected by a bond. Since 
for any dimer covering of the lattice, the bond ij is ei- 
ther occupied by a dimer or not, the number of dimers 
Uij on the bond can only take the values or 1. We 
define p{ij) = {nij) as the average of over all dimer 
coverings. Thus p{ij) is the probability that bond ij is 
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occupied by a dimer. For the star lattice with the same 
dimer weight on all bonds, symmetry dictates that p{ij) 
may take at most two different values as the bond ij is 
varied: p{ij) — pt if ij is a triangle bond and p{ij) = pi 
if ij is a linking bond. Explicit calculations give 

Pt = l, Pi = \ (26) 

Before turning to the derivation of this result, let us show 
that Pt and pi are related to each other. This can be seen 
by considering the average number of dimers in the whole 
system, given by n^j) = Npt + (N/2)pi. Here the 

sum goes over all bonds and the prcfactors of pt and 
Pi are the number of triangle bonds and linking bonds, 
respectively, in a system with N sites. The average is 
trivially equal to N/2 since each dimer covering satisfies 
the close-packed hard-core constraint. Thus 

Pt + lpi = l- (27) 




FIG. 4: The presence of a dimer on bond ij leads to a set 
of "spent" arrows and constraints as shown (the filled circles 
symbolize the constraints). This collection of arrows and con- 
straints is referred to as the dimer's imprint, (a) Imprint of 
a dimer on a triangle bond ij. (b) Imprint of a dimer on a 
linking bond ij. The probability p{ij) of having a dimer on 
the bond ij is found by inserting the number of spent arrows 
and constraints ATVa and AA/'c into Eq. (|3U|I . 



A. Fermionic approach 

The probability p(ij) for having a dimer on bond ij is 
given by^ 

p{ij) = (n,,) - -(^.V,) = ~Grj. (28) 

In this expression the site labels i and j have been chosen 
such that Aij > 0, i.e. the Kasteleyn arrow goes from i 
to j. Eq. (pS)) can be derived3^ from the form of the 
generating function for dimer coverings when the dimer 
weights are allowed to be arbitrary (the dimer generating 
function reduces to our Z when all dimer weights are 
equal). The result now easily follows by using Eqs. 
Q-dlll) to evaluate Eq. (|28|l. 

B. Arrow approach 

The probability p{ij) = (nij) of having a dimer on 
bond ij is given by the ratio 

Pi^j) = ^ii^^piil (29) 

where Z{{nij — 1}) is the number of dimer coverings in 
which bond ij is occupied by a dimer and Z is the total 
number of dimer coverings. In the sections that follow wc 
will discuss other quantities that can be written in this 
way as well, i.e. as the ratio of a numerator, call it Z' in 
general, and the denominator which is Z in all cases. In 
Sec. IIV Al the latter was expressed as Z = 2^"~^<^ where 
Afa is the number of flippable arrows (which are Ising (i.e. 
Z2) degrees of freedom) and Afc is the number of inde- 
pendent constraints on these arrows. For all quantities of 
interest, the numerator Z' can be expressed in terms of 
arrows in a completely analogous way, the only difference 



being that the calculation of Z' corresponds to imposing 
additional requirements on some of these arrows. For 
example, some arrows may be required to point in defi- 
nite directions. If these additional requirements are in- 
compatible with the arrow constraints it follows that the 
numerator vanishes. If, on the other hand, they are com- 
patible with the constraints, they will cause the number 
of remaining constraints and flippable arrows to change 
by AA/'c and AA/'a, respectively, compared to their values 
in the calculation of Z. This gives the result 

£1^2^A^„-AAA.^ (30) 

The calculation of such a ratio therefore amounts to iden- 
tifying the values of AA/'a and AA/'c for the quantity of 
interest. 

Having explained the calculation in general, let us now 
apply this arrow approach to Eq. (f29|) . for which the 
dimer coverings that contribute to the numerator are 
characterized by having a dimer on bond ij. Let us 
first consider the case of ij being a triangle bond. The 
presence of a dimer on such a bond implies that the ar- 
rows on site i and site j must both point into the as- 
sociated triangle. The arrow constraint for the triangle 
then implies that the arrow on the third site of the tri- 
angle must point out of the triangle. Next, as each of 
the three fixed arrows on the triangle also belongs to a 
linking bond, the constraint on each of the three linking 
bonds fixes the direction of the arrow on the other end 
of the linking bond. Thus the directions of six arrows 
become fixed by having a dimer on the triangle bond 
ij. These arrows are therefore no longer flippable, Ising 
degrees of freedom, so AA/'a = —6. Furthermore, four 
constraints are also spent (the triangle constraint and 
three linking bond constraints) so AA/'c = —4. This gives 
p{ij) = 2~^~^~^^ = 1/4 = Pt. This case is illustrated in 
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Fig. Ha). 

Next consider the case of ij being a linking bond. Hav- 
ing a dimer on such a bond implies that the two arrows 
must point into (towards the center of) the linking bond, 
giving AJ\fa = —2. Furthermore, AA/'c = — 1 since the 
constraint on that linking bond is then spent. Thus 
p{ij) = 2~^~(~^) = 1/2 ~ pi. This case is illustrated 
in Fig. lUb). 

The set of spent arrows and constraints that results 
from having a dimer on a given bond will be referred 
to as the dimer's imprint. Fig. IH^a) and (b) shows the 
dimer imprint for a triangle-bond dimer and a linking- 
bond dimer, respectively. Thus from counting the num- 
ber of arrows and constraints in the imprint of a dimer 
on the bond ij, one can use Eq. ([50)1 to calculate the 
probability p{ij) of having a dimer on that bond. 

VII. DIMER CORRELATIONS 

In this section we consider correlations between dimers 
in the star-lattice CDM with equal dimer weights. The 
main focus is on correlations between two dimers, but at 
the end we also briefly discuss correlations between three 
dimers. 

The probability that bonds ij and kl both have a dimer 
is given by p{ij, kl) — {fiijUki). The correlations between 
dimers on bonds ij and kl are quantified by the dimer- 
dimer correlation functioit^ 

c{ij,kl) = p{ij,kl) - p{ij)p{kl) 

= {n-ijnki) ~ (n^j)(nki). (31) 

This quantity can have either sign and vanishes for un- 
correlated dimers. Our results for c(ij, kl) are presented 
in Fig. [5] (all numbers are in units of 1/16). We find that 
a dimer is uncorrelated with all dimers located further 
away than on one of the closest triangles. In particu- 
lar, dimers on linking bonds are uncorrelated with each 
other. 



A. Fermionic approach 

In the fermionic approach {riijUki) is given by a four- 
point correlation function, 

{riijUki) = (ipi-ipj^pk-ipi) . (32) 

This expression for {uijUki) can be derived from the gen- 
erating function of dimer coverings in the same way as 
Eq. ([281) .^'^ (In this derivation the site labels were chosen 
so that the two Kasteleyn arrows go from i to j and from 
k to I, respectively; cf. the remark after Eq. ([^5]) ). Since 
the fermionic action is quadratic, the four-point function 
{ipiipjtpk'4'i) can be evaluated using Wick's theorem. In- 
serting the Wick expansion into Eq. (|3ip and using Eq. 
the term GijGki cancels, giving 

c{ij,kl) = GuG,k-G,kG,i. (33) 




1 2 



FIG. 5: Dimer-dimer correlations on the star lattice. The 
number next to a bond kl is the dimer-dimer correlation func- 
tion c{ij, kl) (in units of 1/16) where the reference bond ij is 
shown in bold, (a) The reference bond is a triangle bond, (b) 
The reference bond is a linking bond. In both (a) and (b) 
only bonds kl with c{ij, kl) ^ are shown. Note that the 
correlations between dimers on different bond types can be 
read off from both (a) and (b). 



The property (|25p of the Green function then implies 
that also the dimer-dimer correlation function will vanish 
beyond a very short distance. (Note that Eq. is not 
valid for the case ij — kl when one instead has {nfj) — 
-Gij, since nf^ = n^j. Thus c{ij,ij) = -Gij - G^-.) 
The results from the Green function calculation of the 
dimer-dimer correlation function are shown in Fig. [5] 

B. Arrow approach 

The probability p(ij, kl) = {uijUki) is given by 

p(.,,feO- ^""'-=^'""^^" (34) 

where Z{{nij = l,nki = 1} is the number of dimer cov- 
erings with dimers on both ij and kl. Eq. (1341) is a ratio 
of the type ((30|) that can be evaluated with the arrow 
representation as explained in Sec. IVIBI One starts by 
laying down the imprints of spent arrows and constraints 
for each of the two dimers (these imprints are shown for 
the two types of dimers in Fig. [3]). If the dimers are 
sufficiently far apart that their respective imprints have 
no arrows in common, the two dimers contribute addi- 
tively to ZWa — AA/'c- This gives p{ij, kl) = p{ij)p{kl), 
i.e. the dimers are uncorrelated.-^ In contrast, dimers 
are correlated if their imprints overlap. In this situation 
p{ij, kl) = if the overlapping parts do not match, which 
happens if the two dimers cannot both be present in a 
dimer covering. Alternatively, if the overlapping parts do 
match, the value of AA/'a — AA/'c is larger (i.e., less nega- 
tive) than the sum of the contributions from each dimer 
considered separately, thus giv'mg p{ij, kl) > p{ij)p{kl). 
This arrow calculation reproduces the results for c{ij, kl) 
shown in Fig. [Sj Note that a dimer on a triangle bond 
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(Fig. [3a)) is correlated with dimers further away than 
a dimer on a hnking bond is (Fig. [Hb)). This difference 
is a consequence of the different sizes of their respective 
imprints shown in Fig. [D In particular, the imprint of 
linking bond dimers is so small that such dimers arc only 
correlated with dimers on the two triangles touching the 
linking bond. Dimers on triangle bonds, having a some- 
what bigger imprint, are correlated only if they are on 
the same or neighboring triangles. 

The notion of an imprint can be generalized to more 
than one dimer. For example, the imprint of two dimers 
at ij and kl is given by the collection of spent arrows 
and constraints resulting from those two dimers. If the 
arrows are mutually compatible, the probability p{ij, kl) 
can be found by inserting the number of spent arrows 
and constraints into Eq. ([50)1 . Another example will be 
given in the next subsection. 



C. Correlations between three dimers 

Correlations between three (or more) dimers can 
be similarly analyzed. Let us briefly consider the 
three-dimer case. This will involve the probability 
p{ij, kl,mn) = {nijUkirimn) that the bonds ij, kl, and 
mn are all occupied by a dimer. Note that even if all 
three dimers are pairwise uncorrelated, i.e. p{ij, kl) = 
p{ij)p{kl), p{ij,mn) = p(ij)p{mn), and p{kl,mn) = 
p{kl)p{mn), it isn't necessarily true that p{ij, kl,mn) 
equals p{ij)p{kl)p(rnn) . As an example, consider the sit- 
uation with the bonds ij, kl, and mn as shown in Fig. 
[SI It follows from the previous analysis in this section 
that dimers on these three bonds are indeed all pairwise 
uncorrelated. However, it is easy to see that if there is a 
dimer on both ij and kl, there must also be a dimer on 
mn. Thus p{ij, kl, mn) = p{ij)p{kl) = 1/16 in this case. 
This result is also found from the fermionic approach, 
using thali^ 

{nijnkinmn) = -{ipt^l^ji'ki'ii'Tni'n) (35) 

which is then evaluated using Wick's theorem (Eq. (1551) 
can be derived in the same way as Eqs. (|28p and (l32p V 
Alternatively, one can use the arrow representation. The 
imprint of spent arrows and constraints for p{ij, kl,mn) 
is shown in Fig. [Sj The fact that this is also the imprint 
for p(ij,kl) shows that p{ij,kl,mn) — p[ij,kl). From 
this figure one sees that AAAa = — 14 and AA/'c = — 10, 
and thus p{ij,kl,mn) = 2-i''-(-io) = 1/16. 



VIII. VISON CORRELATIONS 

The excited-state manifold corresponding to the lowest 
excitation energy of the QDM Hamiltonian ([3|) is spanned 
by states that have eigenvalue —1 for a^{D) on two of 
the dodecagons and -1-1 for all the others. The two lo- 
calized (T^{D) = —1 excitations are called (point) visons. 




FIG. 6: Imprint of spent arrows and constraints for 
p{ij,kl,mn) where the three bonds ij, kl, and mn are as 
shown. This imprint is the same as that for p{ij,kl). It 
has AA/'a ~ —14 and AA^c = —10, giving p{ij,kl,mn) = 
p{ij,kl) = 1/16. 



A vison can thus be labeled by the center of the host- 
ing dodecagon. Let us find an expression for the state 
1^7.,/) (in a given topological sector) with visons at the 
dodecagon sites / and J (we use capital letters for sites 
on the dual lattice; the dodecagon centers form a subset 
of all dual lattice sites). The reasoning is analogous to 
the one used for the kagome lattice Consider a path F 
on the dual lattice that connects sites I and J , as shown 
in Fig. [71 For an arbitrary dimer covering |c), let the 
operator A^r count the number of dimers Ny{c) on the 
bonds intersected by F. Then the operator (—1)''^^ mea- 
sures the associated parity, i.e. whether N-p{c) is even 
or odd. The dimer configuration around a dodecagon D 
corresponds to one of the 32 even-length loops around the 
dodecagon having an alternation of bonds with/without 
a dimer (cf. the discussion around Eq. ^). As (j^{D) 
shifts the dimers by one bond along this loop, it can 
be seen that a^{D) conserves the parity if F crosses the 
loop an even number of times, and flips the parity if F 
crosses the loop an odd number of times. The latter 
case is only realized if D is one of the endpoints / or J 
of F. Thus (7^(i:))(-l)*r (-1)'^^ a"" (D) \i D ^ I,J 
and a'-" {D){-1)^^ = -{-l)^^a''{D) \i D = I, J. It fol- 
lows that the two-vison state we are looking for is given 
by = (-l)^r-|RK) (here |RK) is the Rokhsar- 

Kivelson ground state in the given topological sector) 
since it satisfies (t'=(D)|*/,j) = for D ^ I,J and 

(T^(D)|vI'/ j) = for D = I,J. The vison-vison 

correlation function is equal to (the absolute value of) 
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FIG. 7: Two examples (top and bottom) of a pair of visons on 
the star lattice. The two visons are located at the dodecagons 
whose centers are sites / and J on the dual lattice. A path F 
on the dual lattice (dotted line) connects the two vison sites. 
The vison- vison correlation function is equal to |/c — /o| where 
fe ifo) is the fraction of all dimer coverings that has an even 
(odd) number of dimers on the n bonds intersected by F. The 
top figure is an example of the case considered in Sec. IVIII Al 
with n(= Tlx) = 3. The bottom figure is an example of the 
case considered in Sec. IVIII Bl with = 2, Uy = 1 and 
n — 3. The inset shows the axes and unit cell in Fig. |3] 
translated to this figure. 

the ground-state expectation value of (— 1)^"", 

z;(/,J)^|(RK|(-l)*-|RK)|. (36) 

Clearly v{I,J) — since it is the overlap of |RK) and 
1 4'/.,/) which are eigenstates of the Hamiltonian with dif- 
ferent eigenvalues. As (— 1)^^ is diagonal in the dimer 



basis, v{I, J) can alternatively be expressed as the corre- 
lation function 

«(/,J) = |((-l)^-)| (37) 

in the classical dimer problemi^ i^°i^^ It is instructive to 
rederive the vanishing of v{I, J) from Eq. (|37p. which is 
what we consider in the remainder of this section. 



A. Fermionic approach 

By changing the signs of the matrix elements Aij = 
—Aji for bonds ij intersected by F (which corresponds 
to reversing the direction of the Kasteleyn arrows on 
these bonds), a dimer covering |c) is counted with 
the sign 

(_l)Wr(c) as required by Eq. §7]). In the 
fermionic formulation of this expectation value, these 
sign changes arc implemented by multiplying exp(— S) 
by a compensating factor in the form of the "string" 
nfe-)erexp(2V'»^»jV'i) = nfe)Gr(l + where 
the product runs over all bonds intersected by F. Thus 
the vison-vison correlation function involves the expec- 
tation value of this string, 




We may choose F to only visit dodecagon sites, so that 
all bonds intersected are of the linking bond type. 

As an example of the evaluation of Eq. ([55]) we will 
consider the simple case shown in Fig. Eftop), where 
/ and J have the same y coordinate and are separated 
by a distance n, so that F can be taken as a straight 
line of length n along the x axis. Thus F intersects n 
linking bonds. Let / = {X, Y) and J = {X + n,Y) where 
iX,Y) = {x + l/2,y+l/2). Then v{I,J) = v{X,Y;X + 
n,Y) = \v{n)\ where v{n) is the expectation value of 
the string. A direct evaluation of Eq. (|38p gives v{l) — 
v{2) = 0. For n > 3 we write 



v{n) = / JJ(1 -I- 2ilj^+i^y^eipx+i^y+i^2) ) = v{n - 1) + 2 < 



]^(1 + 2il;x+t^yfiipx+t,y+ia) 



• (39) 



Using Wick's theorem, the last expectation value is the sum of v(n — I) {ipx+n,y,6''l'x+n,y+ 1,2) and 



]^ (1 + 2^px+i^yfi^^:x+l^y+l^2) 



(1 + 2ll}.j;+n-l.,yfi '4^x+n— 1, y +1,2) fpx+n,y,6fpx+n.y +1,2 



(40) 



where the subscript "ext" (external) means that when evaluating this expression using Wick's theorem, the last two 
Grassmann variables (GV's) tpx+n,y,6 and ipx+n,y+i,2 should not be paired with each other. Because of Eq. (f25l) 
any pairings between GV's whose x coordinates differ by 2 or more will vanish. It follows that the contribution to 
(1401) from the first term in 1 -|- 2ipx+n-i,y,6''Px+n-i.y+i.2 will vanish. Furthermore, the only non-vanishing pairings of 



13 



ipx+n,yfi and il^x+n,y+i,2 are with the GV's ')/'x+n-i,y,6 and ipx+n-i.y+1.2 of the neighboring bond. Thus (j40p is equal 

to 2v{n - 2){lljx+n-l,y, l,2^x+n,y,6'4'x+n,y+l,2}cxt^ SO that 

= u(n - 1)[1 + 2{'lp^+n,y,6lpx+n,y+l,2)] 

2>](41) 



The coefScients of v{n — 1) and v{n — 2) in this ex- 
pression vanish, so v{n) = also for n > 3. (In fact, 
we may use Eq. (HTj) to conclude that v{n) = even 
without evaluating these coefficients, as it follows from 
an induction proof, using v{l) — v{2) — 0). Thus 
v{X,Y;X + n,Y) = \i{n)\ = for all n. 



B. Arrow approach 

In this subsection we will use arguments based on the 
arrow representation to consider the vison-vison correla- 
tion function v{X^ Y; X + nx,Y + Uy) with both rix and 
riy arbitrary. In this case F can be taken to be made 
up of two straight line segments of length \nx\ in the x 
direction and length \ny\ in the y direction. An exam- 
ple with Ux and ny positive is shown in Fig. [TKbottom). 
Again, we note that all the bonds intersected by F are of 
the linking bond type. The fraction of dimer coverings 
that has a specific distribution of dimers on the n inter- 
sected bonds (n = \nx\ + \ny\) is equal to 2~". This is 
because the presence or absence of a dimer on a linking 
bond has a probability = 1— Pt = 1/2, and dimers can 
be distributed independently on the linking bonds inter- 
sected by F. This independence can be seen by laying 
down arrows corresponding to a particular dimer distri- 
bution on these linking bonds; the bonds will contribute 
additively to AA/'a — AjVc- (Note that for a linking bond, 
not only the presence, but also the absence of a dimer, is 
associated with a definite configuration of the arrows on 
the two sites connected by that linking bond.) Distribut- 
ing dimers on the intersected bonds in all possible ways 
should generate all dimer coverings. To check that our 
arguments produce this result, we note that the number 
of ways to lay down k dimers on n bonds is (^) and thus 
the total fraction of configurations generated this way is 
Sfe=o (fc) = 1 as expected. Furthermore, the frac- 
tion of configurations fe and /o with an even and odd 
number of intersected dimers, respectively, is given by 

/. - E (:) = \. m 

0<k<n 

- T. ^ \- m 

k odd ^ ' 
0<k<n 

Since the vison-vison correlation function (|37[) can be 
written |/c — /o|, Eqs. imply that it vanishes. 



IX. MONOMER CORRELATIONS 

All properties considered so far have been for a system 
in which every site is touched by exactly one dimer. A 
different but closely related problem involves a lattice in 
which some sites are occupied by monomers, meaning 
that those sites cannot be touched by a dimer. In this 
section we consider a system with two monomers. We 
will be interested in the monomer- monomer correlation 
function, defined as 

m(i,j)^ (44) 

where, for the CDM with equal dimer weights, Zjyiyyi{i^j^ 

is the number of possible dimer coverings when the sys- 
tem has monomers at sites i and j,^^ and Z is the same 
as before, i.e. the number of dimer coverings for the 
system in the absence of monomers. More generally, if 
not all dimer weights are equal, the two quantities in Eq. 
(|44l) are the generating functions for the system with and 
without the monomersi^ 

The monomers in a CDM can be characterized as con- 
fined or deconfined. A deconfined phase is character- 
ized by a nonzero value of 171(1, j) in the limit of infinite 
monomer separation |j — :/| 00, while in a confined 
phase m{i,j) goes to zero in this limit. Using Fisher's 
mapping, ^'^ Moessner and Sondhiii were able to show 
that the star-lattice CDM has a phase transition as a 
function of the dimer weights in which monomers are con- 
fined in one phase and deconfined in the other. They did 
this by relating certain sums of four monomer- monomer 
correlation functions to a spin-spin correlation function 
in the square-lattice Ising model. This approach does 
however not allow one to deduce m{i,j) itself. 

In this section we calculate the monomer-monomer cor- 
relation function for a star-lattice CDM with equal dimer 
weights. We find that m(i,j) — 1/4 for all possible 
choices of the monomer sites i and j (with i ^ j) ex- 
cept when i and j are connected by a linking bond, in 
which case 171(1, j) = 1/2. Thus the monomers are de- 
confined in this case. As we discuss in Sec. IIX C[ these 
results are consistent with the analysis in Ref. illi . 

A. Fermionic approach 

The lattice graph of the system with monomers at sites 
i and j differs from the original lattice graph in that sites 
i and j, as well as all bonds connected to them, are re- 
moved. This modified lattice graph is still planar so the 
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FIG. 8: A particular case of the two-monomer system for which Eq. (|47|l for the monomer-monomer correlation function m{i,j) 
is evaluated in the text. The two monomers are located at sites i and j of the original lattice (shown as open circles). In this 
figure, i = {x,y,6) and j — {x + n,y + 1,2) with n = 6 (cf. the site labeling conventions in Fig. |3]and the inset in Fig. (Tjl. 
The removed bonds connected to these sites are shown as dashed lines. The bonds whose Kasteleyn arrows are reversed in this 
fermionic calculation are shown as thick dotted lines. 



Pfaffian method can again be appliedi^i^ Due to the bond 
removals, the question of signs for the Kasteleyn matrix 
elements needs to be readdressed. Let us consider the 
generic situation when the monomers are sufficiently far 
apart that they sit on separate faces; for the star lat- 
tice each of the faces hosting a monomer has 21 bonds 
along its perimeter. Consider a "path of faces" connect- 
ing the two monomer sites, with all but the two end faces 
on this path being dodecagons. Neighboring faces along 
this path then share one linking bond. It can be seen 
that the original choice of Kasteleyn arrows violates the 
clockwise-odd sign rule around the two end faces host- 
ing the monomers, while the sign rule is respected on the 
other faces on the path. A valid choice of signs is then 
obtained by reversing the Kasteleyn arrows on each of 
the single bonds shared by neighboring faces along the 
path, such that one arrow is reversed around each of the 
two end faces, and two arrows are reversed around each 
of the interior faces. ^'^ An example (which we will return 
to below) is shown in Fig. [S] 

We now give a derivation of the general expression for 
m{i,j) in the fermionic approach (see also Ref. 0). We 
have Zmm{i,j)/Z = |Z„„(i,j)/Z| where, by analogy 
with the original lattice graph, Zmmihj) is given by 



lihj) = y I^V'jj exp(-S'jj). 



(45) 



Here / V^ij ~ / 11; ^V^' ^^^^ Si, 



where A is the modified Kasteleyn matrix (i.e. with the 
necessary sign changes) and the prime on the product 
and sum means that the Grassmann variables for sites i 
and J are excluded. To express Zmmi'i,j) as an expecta- 
tion value with respect to the original lattice graph, we 
use that exp{—Sij) can be written as 

exp ^ ^jfcV'fc^ exp ^ Ajii/ji^ exp(-V'jAijV'j) 

X I n (1 + '2^mA„^n^n) I exp(-5) (46) 
\(mji)eA / 



Here the fourth factor that takes care of the sign rever- 
sals is a "string" of the same type encountered in the 
fermionic calculation of the vison-vison correlation func- 
tion in Sec. IVIII Al the product again goes over the 
bonds whose Kasteleyn arrows are reversed. Further- 
more, / T^'ipij = (up to a possible sign dif- 
ference). The presence of ipiipj here implies that when 
dJS]) is inserted into ([l5|) . the first three factors in ([^5]) 
can be replaced by 1, as the omitted terms will give no 
contribution to Zmmihj) owing to ^lJf = ip'j = 0. Thus^ 



m{i,j) 



\ \(fci)GA / / 



(47) 



We will now evaluate Eq. ([TT]) for the case that the 
left (right) monomer is at site i (j) where i — {x,y,6) or 
{x,y + 1,2) and j = (x + n,y,6) or {x + n,y + 1,2), where 
n > 0. Thus the x coordinates of the two monomers differ 
by n. The directions of the Kasteleyn arrows on the n—1 
vertical linking bonds between i and j are reversed. The 
particular case i — {x,y,6), j — {x + n,y + 1,2) with 
n = 6 is illustrated in Fig. [S] Consider the expression 



Y[{l + 2ljjx+e,y,6 



i'3 



(48) 



whose absolute value is m{i,j). When evaluating this 
using Wick's theorem, all terms involving a 1 from one 
or more of the factors (1 + 2'ijjx+e.y.6'>Px+i,y+i,2) will give 
zero contribution due to the property (pS)) of the Green 
function. Hence the only surviving term in (|48[) is the 
one with all 2n Grassmann variables present. 



Y[ '^4'x+i,y,6''Px+e,y+l,2 



tpj) = rhap{n). (49) 



Here a and (3 are the intra-cell coordinates of sites i and 
j respectively, i.e. a and /? can take the values 2 and 6. 
We find m62(l) = rheeil) = -m22(l) = -m26(l) = 1/4, 
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FIG. 9: Arrow calculation of monomer-monomer correlation function, (a) Imprint (set of spent arrows and constraints) for a 
single monomer located at site i. The open circle denotes the removed arrow on the site. The filled circle denotes the spent 
constraint, (b), (c), and (d): Imprints for some two-monomer configurations. For each of these, m{i,j) can be found by 
inserting the number of spent arrows and constraints into Eq. (|30p . 



while for n > 2, 

rhapin) = 2[mQ,6("- - l)(V'x+n-i,y+i,2V'i) 

= -^[^«,6(" - 1) + ?TJq,2(" - !)]• (50) 

Thus •fhapin) is independent of (3. Defining rhain) = 
TTT-apin), we have 

ma(n) = —rha{n — 1) (f^ > 2), (51) 

with 7716(1) — — TO2(1) = 1/4. It follows that 
the monomer-monomer correlation function m(i,j) — 
\'nT'a{n)\ = 1/4 for any n, including the limit n oo. 

We note that both the monomer-monomer correlation 
function considered here and the vison-vison correlation 
function considered in Sec. I Villi contain a "string" in the 
fermionic formulation. This string causes the number of 
pairings from Wick's theorem to grow exponentially with 
the vison/monomer separation. However, in both cases 
the property (j25p comes to the rescue by making most 
of these pairings vanish. As a result the calculations are 
analytically tractable regardless of the vison/monomer 
separation. This should be contrasted with the situation 
for the triangular-lattice CDM discussed in Ref. [3, for 
which the exponential growth of the number of pairings 
made the fermionic calculation of the monomer- monomer 
correlation function prohibitively difficult for all but the 
smallest values of the monomer- monomer separation. In- 
stead this calculation was there carried out with a deter- 
minant formulation first used by Fisher and Stephenson 
for the square latticei^ This determinant formulation was 
also used for the Pfaffian calculation of the monomer- 
monomer correlation function of the kagomc lattice CDM 
recently reported in Ref. 

Despite the essential simplifications resulting from Eq. 
(j25p . the fermionic approach to the monomer- monomer 
correlation function is still less economical and less intu- 
itive than the arrow approach to be discussed below. For 



some cases of small monomer separation, the fermionic 
calculation of m(i,j) on the star lattice also has an addi- 
tional technical complication: the modified lattice graph 
has cut points, the existence of which violates the condi- 
tions for the explicit proof of the clockwise-odd rule given 
in Ref. |l|. While this complication can be dealt with (the 
Pfaffian method can be used also for a graph with cut 
points^), and one reassuringly finds^^ the same result for 
m{i,j) as obtained from the arrow representation, the 
arrow derivation is much simpler, and therefore we do 
not discuss these fermionic calculations further here. 



B. Arrow approach 

The monomer-monomer correlation function (|44p is yet 
another example of a ratio of the general type discussed 
in Sec. IVIBI that can be calculated using the arrow ap- 
proach. Thus we must again identify the appropriate val- 
ues of AJVa and AJVc in Eq. ([50)1 . Let us start by looking 
at the situation when the two monomer sites i and j are 
far apart. Consider the monomer at one of the sites, say 
i. Its contributions to AA/'a and AA/'c are —2 and — 1, 
respectively. Here, the two spent arrow {Z2) degrees of 
freedom are the (removed) arrow on i and the arrow on 
the site at the opposite end of the (removed) linking bond 
that i was originally one endpoint of; in the presence of 
the monomer this latter arrow must point into the neigh- 
boring triangle. Furthermore, the spent constraint con- 
tributing to AA/'c is the constraint on this (removed) link- 
ing bond.^^ The set of spent arrows and constraints for an 
individual monomer will be referred to as the monomer 
imprint (see Fig. [Hla)). When the two monomers are 
far apart their contributions to AA/'a and AA/'c simply 
add, giving m{i,j) = 2^(^^^(^^)^ = 1/4. This result is of 
course also valid as the separation between the monomers 
is increased to infinity, and is independent of the relative 
orientation between the underlying lattice graph and the 
vector connecting the two monomers. 
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Next, let us analyze what happens when the two 
monomers are moved closer to each other. As long as 
the monomers are not too close, the situation remains 
the same: The two monomer imprints do not overlap, 
nor do they touch the same triangle, so AAfa and AA/'c 
are given by the sum of the separate contributions from 
each monomer. Thus AA/'a — AA/'c is also additive, giv- 
ing m(i,j) = 1/4. Only when the two monomers get so 
close that their relative positions are as shown in Fig. 
injb) do their respective imprints touch the same triangle 
(the middle one in the figure). Now AA/'a and AA/'c are 
no longer additive, but as the additional contributions to 
ZWo and AA/'c are equal, AA/'a — AA/'c is still additive, 
so again m{i,j) — 1/4. This continues to hold when j 
is moved even closer so that i and j are connected by a 
triangle bond, as shown in Fig. M^c). If, on the other 
hand, j approaches i from the side of i's linking bond, 
they will finally end up being connected by this linking 
bond, as shown in Fig. HKd). In this case AAfa = —2 and 
AWc = -1 so m{i,j) = 1/2. 

C. Connection to spin-spin correlations in the 
square-lattice Ising model 

As shown by Moessner and Sondhi^^ (see also Ref. 
113), by using Fisher's mapping^*' between the ferromag- 
netic Ising model on the square lattice and the classical 
dimer model on the star lattice, the spin-spin correlation 
function for the spin model can be related to a sum of 
monomer- monomer correlation functions for the dimer 
model. In our notation, this relation can be written 

{S{x,y)S{x',y')) ^ ^ m{x,y,a;x' ,y' ,a'). (52) 

a, a' —3,4: 

Here S{x,y) is the Ising spin located at position (x, y) 
of the square lattice. In Fisher's mapping, the dimer 
weights on the bonds corresponding to the original square 
lattice (see Fig. El^left)) are 1/ tanh(/3Jij) (where /3 is the 
inverse temperature and Jy is the ferromagnetic Ising ex- 
change parameter) while the dimer weights on the bonds 
internal to the unit cell (see Fig. [Sljright)) are 1. The re- 
lation ([52| implies that the phase transition of the Ising 
model between states with unbroken/broken spin rota- 
tion symmetry maps onto a confinement-deconfinement 
transition in the dimer modelfii The equal-weight dimer 
model considered in this paper is obtained in the zero- 
temperature limit /3 — > cx) of the Ising model, in which 
case the spins are perfectly ordered so the Ihs of Eq. ([5^ 
equals 1. Let us use our results for m{i,j) to evaluate 
the rhs of Eq. ((52|) for this case. When (x,i/) ^ {x\y') 
all four monomer- monomer correlation functions in (j52p 
equal 1/4, so the rhs is 4 • 1/4 = 1. On the other hand, 
when {x,y) — {x',y') the two monomer- monomer cor- 
relation functions with a — a' vanish since with both 
monomers on the same site, no dimer coverings are pos- 
sible (as there is an odd number of sites left to form 
dimers). As the two remaining terms have m(i, j) — 1/2, 



the sum is now 2-0-|-2-l/2 = l. Thus our results for 
m{i,j) do indeed satisfy Eq. ([5^ . 



X. COMPARISON WITH DIMERS ON THE 
KAGOME LATTICE 

As briefly summarized in Sec. IH dimers on the kagome 
lattice have very interesting properties j ^^'^^i^^'^'^i^'^'^^i"^^ 
One may consider the existence of an arrow representa- 
tion of dimer coverings on the kagome lattic o^^i^^'^^ to 
be a "fundamental" property, in the sense that the other 
properties can be derived from it. Alternatively, using 
the Pfafhan method, the "fundamental" property is that 
the Green function of the kagome lattice vanishes beyond 
a very short distancei^i^i 

As shown in this paper, the properties of dimers on the 
star lattice closely parallel those on the kagome lattice. 
Again, the most fundamental properties, from which all 
others can be shown to follow, are the existence of an ar- 
row representation, or, alternatively, the Green function 
satisfying Eq. . (Note that the number of dimer cov- 
erings Z considered in Sec. llVl is also related to the Green 
function: Z ~ Vdet A = l/ydetG.) In the following 
we will discuss the similarities and connections between 
dimers on the star and kagome lattices in some more 
detail. We emphasize that, unless explicitly stated oth- 
erwise, all comparisons in this section between classical 
dimers on these two lattices (including those made above) 
are for the case when all dimers have equal weights, which 
is the case we have considered in this paper. At the end 
we will briefly comment on the star-lattice CDM for un- 
equal dimer weights, whose properties are more conven- 
tional. 

First, however, let us note a close relationship be- 
tween the star and kagome lattices themselves.— i^i^i^ 
These two lattices can be regarded as "expanded" and 
"reduced" versions of each other, in the following sense: 
Starting from the kagome lattice with Nk sites, touch- 
ing triangles can be separated by splitting their common 
site into two sites connected by a new bond. The result- 
ing "expanded" lattice is a star lattice with Ns = 2Nk 
sites (the new bonds are the linking bonds of the star 
lattice). Reversely, starting from a star lattice with Ng 
sites, the linking bonds can be shrunk to zero length and 
the two sites at its ends merged into one. The resulting 
"reduced" lattice is a kagome lattice with half as many 
sites. 

The number of dimer coverings. For a kagome lat- 
tice of Nk sites embedded on a closed surface the num- 
ber of dimer cover ings 13 19121,22^^ Zk{Nk) = 2^^/3+1. 
As shown in Sec. IIV Al the corresponding result for a 
star lattice with sites is Zs{Ns) = 2^»/^+^ Thus 
the kagome and star lattices that are each other's re- 
duced/expanded lattices with Ng = 2Nk have exactly 
the same number of dimer coverings: 

Zg{Ns) ^ ZkiNk). (53) 
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A related special property of the kagome and star lattice 
is that they have very simple expressions for the entropy 
per site/dimer (a simple rational number times log 2), 
in contrast to e.g. the square, honeycomb and trian- 
gular lattice.- The squagomei2i2^ and triangular-kagomc 
latticeSS. also have this property; we will comment more 
on this in Sec. IxTB] 

A consequence of Eq. ((53l) is that the Hilbert spaces 
of QDMs defined on the kagome/star lattices that are 
reduced/expanded versions of each other have the same 
dimension. One can therefore construct one-to-one map- 
pings between the basis states (dimer coverings) in these 
Hilbert spaces. One possible mapping is as follows: Pick 
an arbitrary dimer covering \ck) on the kagome lattice 
and map this to an arbitrary dimer covering |cs) on the 
star lattice. Then map the dimer covering Y\(^^{H)\ck) 
on the kagome lattice to the dimer covering Y\ ^^{D)\cs)- 
Here the product runs over some set of hexagons in the 
kagome case and over the associated set of dodecagons 
in the star lattice case (a given hexagon in the kagome 
lattice is associated with a unique dodecagon in the star 
lattice by the "extension/reduction" procedure described 
earlier). In this way one gets a one-to-one mapping be- 
tween dimer coverings in given topological sectors on the 
two lattices. This mapping is particularly "natural" from 
the point of view of the QDM Hamiltonian Q and its 
kagome analogue considered in Ref. JJi, as they are given 
by a sum of operators on all hexagons/dodecagons. 

Dimer-dimer correlations. On both the kagome and 
star lattice a dimer is uncorrelated with any dimer fur- 
ther away than on a neighboring triangle. For the kagome 
lattice this means that a dimer is correlated with dimers 
that are up to three bonds awayr^iS^iSi while for the star 
lattice, as shown in Sec. IVHl dimers on triangle bonds 
are correlated with dimers up to four bonds away, while 
dimers on linking bonds are correlated with dimers that 
are only up to two bonds away (in fact, only with the 
triangle-bond subset of these, as a linking bond dimer 
is uncorrelated even with dimers on the nearest linking 
bonds that are only two bonds away). Furthermore, com- 
paring Fig. [5fa) with the corresponding figure for the 
kagome lattice (the latter can be constructed from the 
results in Ref. 24, or calculated from the arrow repre- 
sentation), one finds that the numbers on the triangle in 
the middle are identical while those on the three outer 
triangles have the same magnitude but opposite sign (the 
linking bonds are of course not there in the kagome lat- 
tice figure). 

Vison-vison correlations. The vison-vison correlation 
function v{I, J) vanishes for both lattices (here / and 
J are hexagon/dodecagon sites, respectively). In both 
cases this is most easily seen from the fact that, as dis- 
cussed in Sec. IVIIIl the state with visons at / and J is 
an eigenstate of the QDM Hamiltonian, and the absolute 
value of the overlap between this two-vison excited state 
and the ground state is just the vison-vison correlation 
function. 

Monomer-monomer correlations. The monomer- 



monomer correlation function m{i,j) on the kagome lat- 
tice equals 1/4 for all monomer sites i and j.2P'.'^^ As 
shown in Sec. IIXI the same result holds for the star lat- 
tice, with one exception: when i and j are connected by 
a linking bond, m{i,j) = 1/2. 

Finally, we note that the result (P5|) for the star-lattice 
Green function is limited to the case of equal dimer 
weights that we have considered in this paper. As a con- 
sequence, the very special properties of the star-lattice 
CDM for equal dimer weights do not carry over to un- 
equal dimer weights. In contrast, the special properties 
of dimers on the kagome lattice are more robust, as the 
result (pS)) continues to hold also if the dimer weight is 
allowed to depend on the orientation of the bondi^i^i 



XI. GENERAL FISHER LATTICES AND THEIR 
"REDUCED" LATTICES 

It is interesting to ask whether the results obtained 
here for the star lattice can be generalized to other lat- 
tices. A natural class of lattices to consider in this re- 
spect are those hosting the dimers in Fisher's mapping 
from an Ising model to a dimer problem, of which the 
star lattice is one particular example. Any such "termi- 
nal" lattice in Fisher's mapping will here be referred to 
as a general Fisher lattice, or just Fisher lattice for short. 
The sites (vertices) in a general Fisher lattice can have 
coordination number (henceforth called degree) 1, 2, or 
3. Sites of degree 3 are part of a triangle of three sites 
connected by bonds which we call triangle bonds. All 
other bonds in the Fisher lattice will be called linking 
bonds. 



A. Arrow representation and the number of dimer 
coverings for general Fisher lattices 

An arrow representation of dimer coverings exists also 
for a general Fisher lattice. Again, the arrows are Ising 
degrees of freedom living on the lattice sites, and each 
triangle and each linking bond contribute an arrow con- 
straint. Using the arrow representation, one can calculate 
various properties of dimers on general Fisher lattices, 
analogously to what we have done for the star lattice. 
As an example, we will show that the arrow represen- 
tation gives the correct answer for the number of dimer 
coverings on an arbitrary Fisher lattice. We do this by 
first using Fisher's mapping to calculate this number, 
and then we show that the arrow representation gives 
the same result. 

Thus let us consider Eq. (11) in Ref. [T^, which relates 
the Ising model partition function (defined in Fisher's 
Eq. (3)) to the generating function of dimer coverings 
on the associated Fisher lattice. We are interested in 
the case of equal dimer weights, i.e. Vij = 1, for which 
this generating function reduces to the number of dimer 
coverings Z on the Fisher lattice. This gives (in our nota- 
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tion) Zj = 2^'~^'-''Z. Here Zj is the partition function 
of the Ising model and Nj and Nj^i, are the number of 
sites and bonds in the lattice on which this Ising model 
is defined (which we will refer to as the Ising lattice) . As 
the equal- weight case maps to the zero-temperature limit 
of the Ising model, Zj — 2 since only the two perfectly 
ordered states related by time reversal contribute to Zj 
in this limit. Thus 

Z = 2^'-^-^'+\ (54) 

Next, let us calculate Z using the arrow representation. 
As in Sec. HV^ we have Z = 2^»-^'= where A^^ is the 
number of Z2 arrow degrees of freedom and Afc is the 
number of independent arrow constraints. For a general 
Fisher lattice we have A/'a = Np — N^j^'' where Np is 

the number of sites on the Fisher lattice and Np'' is the 
number of sites of degree 1. The reason for subtracting 
Np^ is that an arrow on a site of degree 1 has to point 
towards the single bond connected to the site (because 
the site has to be touched by a dimer) and is therefore 
not a Z2 degree of freedom. The number of independent 
arrow constraints is Afc = Nt + Ni — 1 where Nt is the 
number of triangles and Ni is the number of linking bonds 
in the Fisher lattice. As in Sec. IIV Al the term —1 in Afc 
comes about because one of the Nt + Ni constraints can 
be deduced from the others. We show in Appendix IB II 
that 

{NF-N'p'^)-{NT + Ni)^Ni^b~Ni (55) 

which means that the result ([54]) is indeed reproduced by 
the arrow representation. In Appendix IB II we also show 
that 

Ni^h-Ni^^{NF-NP)-^N^p'\ (56) 
For a general Fisher lattice we therefore have 

For Fisher lattices with no sites of degree 1 or 2 this 
simplifies to 

Z = 2^^/6+1. (58) 

This simplification applies e.g. to the star lattice graph 
considered in Sec. IIV Al as all its sites have degree 3; 
the result for Z is indeed the same as found there. In 
Appendix IB 21 we check the result (|57|) for some Fisher 
lattices with sites of degree 1 and/or 2. 

B. The "reduced" lattice of a Fisher lattice 

The kagome lattice is the "reduced" lattice of the star 
lattice in the sense described in Sec. [X] An analogous 
reduction procedure, consisting of shrinking the linking 
bonds to zero and merging all sites that come together 



in this process, can be carried out for any Fisher lattice 
F. Unless F is trivial by not containing any triangles, the 
sites in the reduced lattice R so obtained have coordina- 
tion number four. 

As a first example, consider the Fisher lattice asso- 
ciated with the triangular-lattice Ising model, which can 
be shown to consist of triangles, octagons, and "16-gons." 
The reduced Fisher lattice R is a squagome latticej^i^ 
This is a lattice of corner-sharing triangles and therefore 
has an arrow representation of its own as noted in Refs. 

[Ulli. 

As a second example, consider the kagome-lattice Ising 
model. The lattice E (this notation is explained in Ap- 
pendix |B]) is then the star lattice. The Fisher lattice F 
is obtained by replacing each site of the star lattice with 
a triplet of sites connected by a triangle of additional 
bonds. The reduced lattice R of this Fisher lattice is 
the triangular-kagome lattice recently considered in Ref. 
'26'. Although this is not a lattice of corner-sharing tri- 
angles (as its triangles also share edges), it still has an 
arrow representation. To see this, note that any dimer 
on this lattice can be uniquely associated with an "outer 
small" triangle. (Each "big" triangle consists of 4 small 
triangles: 3 outer and 1 inner). Therefore an arrow repre- 
sentation exists in which each arrow must point into one 
of the two outer small triangles that its site connects, 
i.e. the arrows are again Z2 variables. The fact that 
the inner small triangles have a different status than the 
outer small ones is refiected in the fact that the former 
triangles are not present in the Fisher lattice F but only 
appear after F has been reduced to R. 

The number of dimer coverings on F and R are closely 
related, as we have already seen for the star and kagome 
lattice in Sec. |Xl The reduction of F to R involves remov- 
ing one site for each linking bond, and therefore the dif- 
ference Afa—Afc is the same in F and R, and thus so is the 
number of dimer coverings (cf. Eq. ([53])). Since the num- 
ber of sites on the two lattices are related by Np = 27V/j, 
the entropy per site on R in the thermodynamic limit is, 
from Eq. ([55)l . given by sj^ = (1/3) log 2, which agrees 
with previous calculations in the literature i"'^^i^^'^^i^^''^^i^^ 
The fact that the entropy per site is the same for the 
kagome, squagome, and triangular kagome lattice can 
thus be understood to be a consequence of the fact that 
each of these lattices is the reduced lattice of a Fisher 
lattice for which Eq. ([55]) holds. 

Based on the existence of the arrow representation as 
well as recent studies of what we have here dubbed re- 
duced Fisher latticeS j^^i^^'^^ we expect that also other 
properties of dimers on general Fisher lattices and their 
reduced lattices should in general be closely related and 
resemble those of the star and kagome lattice (one excep- 
tion is the purely one-dimensional Fisher lattices consid- 
ered in Sec. [B2|) . 
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FIG. 10: In this figure filled circles represent lattice sites. Left: 6-site unit cell for the star lattice. An Ising pseudospin a — ±1 
is associated with the vertical linking bond connecting the two triangles. Middle: The unit cell graph that results from deleting 
the bonds that, for the given value of a, cannot be occupied by a dimer. (The thick line in the a = +1 graph is the dimer on 
the linking bond.) Right: The decoupled horizontal loops resulting from deleting the appropriate bonds in all unit cells for a 
particular pseudo-spin configuration (the pseudo-spins are shown as signs ±) in a system with L — 4. (For ease of drawing, 
the vertical dimers in the a = +1 unit cells are not shown and the bonds connecting different unit cells have all been made 
horizontal.) 



XII. FINAL REMARKS 

We conclude by mentioning another lattice recently 
considered in the literature that also has an arrow repre- 
sentation. This lattice, which is a hybrid of the kagome 
and star lattices, is shown in Fig. 1 (right) in Ref. fill , 
which considered a dimer problem on this lattice. Using 
the arrow representation one can e.g. reproduce the re- 
sult for equal- weight dimers (u = 1) obtained in Ref. [ll| 
with the Pfaffian method. 
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APPENDIX A: AN ALTERNATIVE 
PSEUDO-SPIN REPRESENTATION OF DIMER 
COVERINGS ON THE STAR LATTICE 

Recently Dhar and Chandra introduced an alternative 
pseudo-spin representation of dimer coverings on lattices 
of corner-sharing triangles.'''^ Wu and Wang used this rep- 
resentation to give a derivation of the number of dimer 
coverings on a finite kagome lattice.— In this Appendix 
we generalize this pseudo-spin representation to the star 
lattice and use it to give an alternative simple derivation 
of the result (jlSp for the number of dimer coverings on a 



finite star lattice on a torus. 

Consider the star lattice as shown in Fig. [TJ Here we 
will choose the 6-site unit cells to be oriented as shown 
in Fig. [TU] (left panel). A unit cell will be labeled (kj) 
where k = 1, . . . , K and £ = 1, . . . , L. Here L is the num- 
ber of horizontal rows of these unit cells and K is the 
number of unit cells in each such row. We will consider 
the star lattice graph to have a torus geometry by con- 
necting opposite ends in the horizontal and vertical di- 
rection. The number of unit cells is where N — 6KL 
is the total number of sites in the lattice. 

We now associate an Ising pseudo-spin variable ake = 
±1 with the vertical linking bond in unit cell {k,£). We 
define aki — +1 (—1) if this linking bond is occupied 
(not occupied) by a dimer. Let us consider how many 
dimer coverings are associated with a given pseudo-spin 
configuration {aki}- One sees that for a given value of 
ake some bonds in the unit cell {k,£) cannot be occu- 
pied by a dimer. By removing these bonds in all the 
unit cells one obtains a new lattice graph that consists 
of L disconnected horizontal closed loops of bonds. This 
is illustrated in Figs. [TO] (middle and right panels). If 
all of these disconnected loops contain an even number 
of sites, the number of dimer coverings associated with 
the pseudo-spin configuration is 2^, since there are two 
dimer coverings for each loop. On the other hand, if some 
loops have an odd number of sites, no dimer coverings are 
possible for the given pseudo-spin configuration. 

It remains to find the number of pseudo-spin configu- 
rations that are consistent with dimer coverings. To this 
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end, we define the variabless* 

K 

Ti^WaM, e=l,...,L. (Al) 
fc=i 

If (Tkg = +1 the unit cell {k, £) contributes 2 sites to each 
of the two loops above and below it, while if aki = ~1 
the unit cell contributes 3 sites to each of these loops. In 
order for the loop going between unit cell rows i and i+l 
to contain an even number of sites, an even number of the 
unit cells contributing sites to it must have pseudo-spin 
(T = — 1. This constraint can be written 

TeTe+i ^ 1. (A2) 

There are L such constraints, one for each loop. However, 
only L — 1 of these are independent, since due to cr^^ = 1 
we have 

{riT2){T2T3) . . . (tlt^) ^ 1 (A3) 

so one of the constraints can be deduced from the oth- 
ers. The number of pseudo-spin configurations that are 
consistent with dimer coverings is therefore 2^/6-(i-i). 
Since each of these configurations is associated with 2^ 
dimer coverings, the total number of dimer coverings is 

Z = 2^/6-(i-i) . 2^ = 2^/*^+! (A4) 

in agreement with Eq. psp . 

APPENDIX B: DETAILS FOR SEC.IXll 
1. Derivation of Eqs. (f55|) - (f56|) 



Next, Fisher introduces a "terminal" lattice, which we 
will call the Fisher lattice F. It is obtained from E by 
leaving sites of degree 1 as they are, replacing each site 
of degree 2 by a pair of new sites (connected by an ad- 
ditional bond), and replacing each site of degree 3 by a 
triplet of new sites (connected by a triangle of additional 
bonds). The number of sites in F is therefore 

Nf = N^E^ + 2N^E^ + 3iVf ' = 2Ni^b + ^{2q- 6)7vj'^ 

(B5) 

Furthermore, the number of bonds in F is 

NF,b = NE,b + n'^e^ + iNf (B6) 

since each site of degree 2 in E gives one new bond and 
each site of degree 3 in E gives 3 new bonds (the triangle 
bonds). The number of triangles in F is 

iVT = 7V|^^ =ivf^-t-^(g-2)7vj«) (B7) 

q>i 

The number of linking bonds in F is the sum of the first 
two terms in NE.b, i-c 

Ni ^NE,b + N^E^ = iV,,b + ^ (9 - 3)ivj«^ + 7vj2) . (B8) 

q>4 

Eqs. ([55]) and (|56p now follow from these expressions and 
the fact that = Tvj^^ and = 2N'^^\ 



Consider an Ising lattice I composed of Nj' sites of de- 
gree q, where q runs over the positive integers. The total 
number of sites and bonds in this lattice are, respectively, 

Nj = (Bl) 

Fisher^*' introduces an "expanded" lattice E in which 
each vertex of degree q > 4 in I is replaced by a "cee" 
of g — 2 vertices of degree 3 and g — 3 additional bonds 
connecting these. Thus the number of sites and bonds in 
E are, respectively, 

Ne - ^7Vl^)+5](g-2)Arj^\ (B3) 

9=1 q>i 

NE,b = Af/,, + 5](g-3)iVj«' = i^g7V(?\(B4) 

where in the last expression Np' = Nj, Np' = Nj , 
and nI^^ = TVj^^ -I- Eg>4(9 " 2)7Vj'^ (by construction, 
N^^^ = for 9 > 4). 



2. Check of Eq. (|57[) for Fisher lattices with sites 
of degree 1 and/or 2 

As a check of Eq. (|57p when Np^ and/or Np^ are 
nonzero, let us consider an Ising model in one dimension. 
First consider the case of periodic boundary conditions. 
All sites in the Ising lattice I are then of degree 2. The 
expanded lattice E is therefore identical to I. Continuing 
to the Fisher lattice F, each site in E is replaced by 2 sites 
in F. Thus F is simply a ID chain with an even number 
of sites with periodic boundary conditions. Obviously, 
2 dimer coverings are possible on F in this case, and 
this indeed also follows from Eq. (|57p when we insert 

Nf = N^F^ and N^^^ = 0. 

Next, consider a ID Ising lattice with open boundary 
conditions. Now the two end sites in I have degree 1 while 
all others have degree 2 as before. Again E is identical to 
I, while F has twice as many degree-2 sites (so the total 
number of sites in F is even). Obviously, only 1 dimer 
covering is possible on F in this case, and this is also 
what Eq. (|57|) gives upon inserting Nf — Np^ + Np'^ 
and N^^'^ = 2. 
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